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Abstract 

In theories with gauge mediated supersymmetry breaking, the scalar tau, 

(fi) is the lightest observable superpartner for part of the parameter space. 

At LEP 2, the production of such a fi pair and their subsequent decays give 

rise to a pair of r leptons plus missing energy from the unobserved gravitinos. 

The angular distributions of the r's are different from those arising from the 

production and decay of W pairs, and thus will constitute an interesting signal 

for supersymmetry. We also consider fi pair production in the complementary 

part of parameter space where the lightest neutralino is lighter than the fi . 
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In spite of many interesting features of supersymmetric theories, the mechanism of super- 
symmetry breaking and how it is communicated to the observable sector has been a major 
area of concern for over a decade. In most of the previous work it has been assumed that 
the supersymmetry is broken in a hidden sector at a scale of ~ 10^^ GeV, and communi- 
cated to the observable sector via the gravitational interaction. However, this scenario has 
a major problem involving flavor changing neutral current (FCNC). This can be avoided if 
supersymmetry is broken in a hidden sector at a scale 10^ GeV, and communicated to the 
observable sector by means of the Standard Model (SM) gauge interaction The number 
of parameters is also smaller in these gauge mediated models (GMSB) which makes them 
more predictive. These models have a very distinctive feature, the gravitino is the light- 
est supersymmetric particle(LSP) and all super particles must ultimately decay to it. This 
feature gives rise to some new signals in colliders that may be observed in the near future. 

One interesting aspect of gauge mediated models is that the role of the next-to-lightest 
supersymmetric particle (NLSP) can be taken by either the lighter scalar tau (fi) or the light- 
est neutralino. In the past year a lot of work has been done |0-|TT|]in these models motivated 
in part by the fact that, when the neutralino is the NLSP, the CDF event (ee77-|-missing 
energy can be explained. However there is a vast region of parameter space where the 
lighter stau can be the NLSP and the lightest neutralino is the NNLSP(Next to Next to 
LSP). Recently, we considered neutralino pair production in these scenarios where the stau 
is the NLSP and showed that the neutralino pair gives rise to 4r + missing energy in the 



final state without any standard model background [T^. In this paper we look at direct 
production of stau in the scenarios where the lighter stau is the NLSP as well as in the sce- 
narios where neutralino is the NLSP. In the first case, since the only decay mode available 
to a lighter stau is to decay into a r and a gravitino, the signal is 2 r + missing energy. 
We compare the signal with the background generated from W pair production (each W 
can decay into a r and a tau neutrino(z/T-)), and show that the angular distribution of the 
individual r's will make the signal look very different from the background. For the 2nd 
case the signal is 2 r +27 with missing energy since each stau decays into a neutralino and 



a r and the neutralino then decays into a hard photon and a gravitino. So the signal itself 
can discern the two cases. We also discuss the angular distribution and energy distribution 
of the decay products in the 2nd case. 

In GMSB models the superparticle masses depend on five parameters M, A, n, tan f3, fi. 
M is the messenger scale, M = A < s >, where < s > is the VEV of the scalar component 
of the hidden sector superfields, and A is the Yukawa coupling. The parameter A is equal to 
< Fs > I < s >, where < > is the VEV of the auxiliary component of s. Fg can be ~ F 



1^ , where F is the intrinsic SUSY breaking scale. In GMSB models, A is taken around 100 
TeV, so that the colored superpartners have masses around a TeV or less. The parameter n 
is fixed by the choice for the messenger sector. The messenger sector representations should 
be vector like (for example, 5 + 5 of SU{5), 10 + 10 of SU{5) or 16 + 16 of SO(IO)) so that 
their masses are well above the electroweak scale. They are also chosen to transform as a 
GUT multiplet in order not to affect the gauge coupling unification in MSSM. This restricts 
n(5 + 5) < 4, n(10 + 10) < 1 in SU{5), and n(16 + 16) < 1 in 50(10) GUT for the messenger 
sector (one 10 + 10 pair corresponds to n{5 + 5)=3). The parameter tan/3 is the usual ratio 
of the up (Hu) and down {H^) type Higgs VEVs. The parameter /i is the coefficient in the 
bilinear term, /iHuHd in the superpotential, while another parameter B is defined to be the 
coefficient in the bilinear term, BfiHuH^ in the potential. In general, fi and B depend on 
the details of the SUSY breaking in the hidden sector. We demand that the electroweak 
symmtery is broken radiatively which determines /i^ and B in terms of the other parameters 
of the theory. Thus we are left with five independent parameters, M, A, n, tan f3 and sign(/i). 
The soft SUSY breaking gaugino and the scalar masses at the messenger scale M are given 



by 0,0 



and 



faHM) = 2 („) / (A) Y. k, C. (2) 



where i = 1,2,3 are the three SM gauge couphngs and ki = 1, 1,3/5 for SU(3), SU(2), 
and U(l), respectively. The Ci are zero for gauge singlets, and 4/3, 3/4, and (V/2)^ for the 
fundamental representations of SU (3) and SU (2) and f/(l)y respectively (with Y defined by 
Q = I3 + Y/2). Here n corresponds to n{5 + 5). g{x) and f{x) are messenger scale threshold 
functions with x = A/M. 

We have calculated the SUSY mass spectrum using the appropriate RGE equations 
T6| with the boundary conditions given by equation (1) and (2), and varying the five free 



parameters. Although in principle the messenger scale is arbitrary (with M/A > 1), in our 
analysis we have restricted 1 < M/A < 10^ and chosen A ~ 100 TeV. For the messenger 
sector, we choose 5 + 5 of SU(5), and varied n{5 + 5) from 1 to 4. In addition to the current 
experimental bounds on the superpartner masses, the rate for 6 — S7 decay puts useful 
constraints on the parameter space [|6|,[7|,p!7| . It is found that for n = 1 and low values 
tan/? (tan/5 < 25), the lightest neutralino xo is the NLSP for M/A > 1. As tan/5 increases, 
fi becomes the NLSP for most of the parameter space with lower values of A. For n > 2, 
fi is the NLSP even for the low values of tan/? (for example, tan/? ^ 2), and for n > 3, fi 
is again naturally the NLSP for most of the parameter space. In Tables 1 and 2 we show 
five sets of spectrum where the lighter stau is the NLSP and five sets where the lightest 
neutralino is the NLSP. We use these scenarios for detailed calculations. 

Let us first discuss the production cross section of the fi pair in the scenarios 1-5 of 
Table 1, where the lighter stau is the NLSP. The total cross-sections are given in Table 3 for 
three LEP2 energies, y/s =172, 182 and 194 GeV . Each of the produced fi will decay into 
a T and a gravitino with essentially a 100% branching ratio. Thus, from fi pair production, 
we obtain final states with two r 's and missing energy. For example, the number of events 
in scenario 5 for y/s =182 GeV, with 100pb~^ luminosity, is 27, and for ^/s =194 GeV , with 
250pb~^ luminosity is 70. The decay, fi tG is fast enough so that it takes place inside 
the detector. (If is much larger than a few 1000 TeV [|18[ , then fi will decay outside 
the detector. In that case the signal will be 2 heavy charged particles passing through the 
detector). 



There is a considerable background from W production, where each W decays into a r 
and a neutrino. However since the staus are scalar particles, their angular distribution is 
different from the W distribution. This shows up in the individual r angular distributions. In 
Figure 1 we show these angular distributions of the rs. As expected the angular distribution 
of the T coming from stau decay is basically flat; the same distribution from W decay has 
signiflcant 6 dependence (as measured from the beam direction). Consequently, the signal is 
almost twice the background for t~ when cos^ is negative. The opposite thing happens for 
r+ i.e the positive values of the cos6' regions need to be searched for r~*" in order to extract 
the signal. 

The information in Fig.l, and its mirror image for the r of the other charge, is made 
more precise in Table 4 where we give the fraction of events that should be found in each 
bin of width 0.2 in cos^+ and in cos6^. 

Figure 2 shows the distribution in the angle between the two tau is very similar in stau 
or W production. 

In Table 3 we also show the cross-section for charged Higgs production in a usual SUSY 
model. The charged Higgs in the GMSB model are too heavy to be produced at LEP 2 but 
in other models they could be light and, because they decay into a r + z/^ with essentially 
100% branching ratio if tan/5 is greater than 2, and because they are produced with the 
same angular distribution as the stau, they are indistinguishable from stau production of 
the same mass in GMSB except through the total production rate which is a little larger than 
for stau. For example the number of events, when the charged Higgs mass is comparable 
to the lighter stau mass in the scenario 5, is 34 for y/s =182 GeV with 100pb~^ luminosity 
and is 89 for ^/s =194 GeV with 250pb^^ luminosity. So if a two tau signal is observed at 
LEP 2, with a limited number of events so that the difference in the production rate is not 
sufficient to discern between these two models, then other signals, such as 4 r + missing 
energy from neutralino pair production |T^, will have to be used. We note here that in 
the case of charged Higgs pair productions, with subsequent decay into r will have similar 
angular distribution as those coming from stau. Thus, the signal for the charged Higgs can 
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also be distinguished from the W-background using the distribution shown in Figl. 

We now discuss the scenarios of Table 2 where the lightest neutralino is the NLSP. In 
Table 5 we give the cross section where staus are produced. The signal in this case is different 
than the signal in the case where stau is the NLSP; here each stau will decay into a r and 
a neutralino and the neutralino will then decay into a gravitino and a photon giving rise to 
2 r + 2 7 + missing energy in the final state. The photons are hard and they can be easily 
seperated from bremsstrahlung. The number of events in scenario 6 for y/s —182 GeV, with 
lOOpb'^ luminosity, is 30 and for ^/s =194 GeV , with 250pb'^ luminosity is 77. Since the 
neutralino is the NLSP in these scenarios, they can be pair produced more easily than the 
staus. Consequently the signal of 2 photons(hard) and missing energy in the final state 
coming from the neutrahno decays could be used to discover SUSY and a detailed analysis 
of the stau production will then give information on particular models. 

We now give three graphs showing the angles and energies of the final state rs and 
photons. There is nothing very surprising in these plots but we include these because they 
might be useful in verifying that detected rs and photons are from this process. In Figure 3 
we plot the angular distributions of the one of the emitted photons and the distribution of 
the angle between the two emitted photons. The angular distribution of the single photon 
is almost isotropic while the angular distribution of the angle between the photons has 9 
dependence; the distribution is large when cos^ is negative. The angular distribution of one 
of the T, or of the angle between the two r looks very similar to Fig. 3. In Fig. 4 we show the 
angle between a tau and a photon from the same stau decay and between a tau and a photon 
which come from opposite stau decays. In Fig. 5, we also plot the energy distributions of 
the decay products, ie., the total missing energy, the energy of one of the photons and the 
energy of one of the r's. Since the massive neutralino carries more than half of the energy 
of the lighter stau, each gravitino carries more than 1/6 of the total energy and the total 
missing energy distribution becomes a maximum at more than 1/3 of the beam energy. The 
energy distribution for the r's and photons maximize at much smaller energy. 

In this paper we concentrated on the production of staus. There are other scalar particles, 
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selectron and smuon, whose lighter component's masses can be less than the beam energy 
but this is usually not the case for LEP II. Even when they are lighter their masses are 
greater than that of the stau and close to the beam energy so that their production is 
suppressed relative to that of the stau. Production of selectron is further suppressed by the 
additional contribution of neutralino exchange in the t-channel while the production rate 
of the smuons is the same except for the effects of the different (larger) mass. There are 
three possible signals for the production of selectrons (smuons): 1) Eb > nie > m^o > rrif, 
where Eb is the beam energy, then the e will decay to e + x° followed by x° — > r + f and 
f — >■ tG. The final state is e+ + + 4r plus missing energy. 2) If Ei, > rrie and m^o > mg 
then the selectron decays directly to an electron and a gravitino; the final state is just e^e~ 
plus missing energy. 3) If Eb > rrii > ruf > m^o then the selectron will decay to e + x° with 
X°— > 7 + G; the final state is e+e^ + 27 plus missing energy. 

In conclusion, we have discussed the production of scalar staus and their decay modes 
in gauge mediated models in the scenarios where the lighter stau is the NLSP as well as in 
the scenarios where the neutralino is the NLSP. We have discussed the background for the 
individual cases and have shown how the signal can be extracted from the SM background. 
We also have noted that the charged Higgs in the supergravity models can have signatures 
that are almost identical to that of the stau in these models so that observation of this signal 
is a sign of new physics but not necessarily of gauge mediated supersymmetry breaking. 

After finishing this work, we came across a preprint by S. Ambrosanio, G. D. Kribs and 



S. P. Martin, |hep-ph/ 970321 1| , which has partial overlap with our work. 
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Texas at Austin. He wishes to thank Duane Dicus of the Center for Particle Physics for very 
warm hospitality and support during his stay there. This work was supported in part by 
the US Department of Energy Grants No. DE-FG013-93ER40757, DE-FG02-94ER40852, 
and DE-FG03-96ER-40969. 
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TABLE CAPTIONS 



Table 1 : Mass spectrum for the superpartners in some scenarios where the lighter stau is the 
NLSP. Note that the 1st and 2nd generation superpartner masses are almost same. 

Table 2 : Mass spectrum for the superpartners in some scenarios where the lightest neutralino 
is the NLSP. 

Table 3 : For each scenario in Table 1 and three beam energies the first line gives the total cross 
section for stau pair production, while the 2nd line gives the total cross section for 
charged Higgs production. 

Table 4 : The double angular distribution of the two tau. The numbers are the per cent of 
total events where each tau is in a bin in cos^ of width 0.2. The first number in 
each bin is for r from e^e^ — > fifi — > ttGG] the number in brackets is for r from 
e'^e" — > WW — > ttUt-Vt. The table is symmetric about the complementary diagonal, 
ie., the numbers for (cos^+ , cos^_) are the same as those for (- cos^_ , - cos^+). 

Table 5 : For each scenario in Table 2 and three beam energies the total cross section for stau 
pair production is shown. 

FIGURE CAPTIONS 

Fig. 1 : The angular distribution of r~, relative to the electron beam axis, from the stau decay 
(dashed line) and from W decay (solid line) . The distribution of the r"*" is the mirror 
image of this graph. 

Fig. 2 : The relative number of events as a function of the angle between the t+ and the 
T~ . The dashed hne corresponds to the rs from the stau decay and the sohd hne 
corresponds to the rs from the W decay. 

Fig. 3 : The angular distribution of one of the photons (dashed line) and the angle between 
the two photons (solid line) from neutralino decay, when the neutralino is the NLSP. 
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Fig. 4 : The angular distribution between a r and a photon. The sohd hne corresponds to the 
T and the photon from the same stau and the dashed hne corresponds to the same 
particles from opposite stau. 

Fig. 5 : The distribution of the missing energy (solid line), the energy of one of the photons 
(dashed), and the energy of one of the r (dot-dashed), when the neutralino is the 
NLSP. 
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Table 1 





Scenario 1 


Scenario 2 


Scenario 3 


Scenario 4 


Scenario 5 


masses 


A = 63.7 TeV, 


A = 33 TeV, 


A = 60 TeV, 


A = 59.7 TeV, 


A = 28 TeV, 


(GeV) 


n=l, M = 4A 


n=2, M = 20A 


n=l, M = lOA 


n=l, M = lOA 


n=2, M = 40A 




tan /3=31.5 


tan/3=20 


tan /9=31.5 


tan/3=28.5 


tan/?=18 




121 


117 


120 


120 


114 




366 


318 


356 


364 


278 


mA 


357 


308 


347 


355 


266 


m^o 


85 


87 


80 


80 


72 


m^i 


158 


156 


149 


148 


128 


m^2 


350 


286 


345 


343 


249 


m^3 


364 


309 


358 


356 


275 


m^± 


157,367 


155,312 


149,361 


127,277 


158,367 




74,249 


73,192 


65,240 


74,236 


65,167 




120,236 


96,184 


116,225 


115,224 


85,159 




664,727 


515,588 


607,673 


605,672 


432,505 




698,740 


558,586 


641,686 


643,683 


472,497 




737,765 


580,601 


683,710 


679,709 


490,508 


ai,2 


735,769 


580,606 


681,715 


678,711 


490,514 


nig 


565 


587 


533 


530 


498 




-343 


-278 


-337 


-336 


-240 
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Table 2 





Scenario 6 


Scenario 7 


Scenario 8 


Scenario 9 


Scenario 10 


masses 


A = 42 TeV, 


A = 44 TeV, 


A = 41 TeV, 


A = 38 TeV, 


A = 42 TeV, 


(GeV) 


n=l, M = 60A 


n=l, M = 80A 


n=l, M = lOOOA 


n=l, M = 800A 


n=l, M = 1200A 




tan/3=21 


tan/3=18 


tan 


tan /3=12 


tan /3=8 




112 


113 


111 


109 


109 


niifi 


280 


300 


296 


277 


316 


rriA 


268 


290 


285 


265 


306 


m^o 


54 


57 


53 


48 


53 


m^i 


98 


104 


96 


87 


98 


m^2 


255 


270 


263 


244 


274 


m^3 


272 


286 


281 


263 


293 


m^± 


98,275 


290,103 


284,96 


266,86 


96,295 




62,174 


73,180 


78,171 


77,157 


90,171 




89,164 


92,172 


92,165 


86,154 


94,170 




419,497 


432,515 


377,472 


446,355 


380,483 




442,473 


484,512 


416,442 


391,412 


448,432 




466,485 


485,506 


436,457 


406,425 


445,467 


nij 

ai,2 


465,492 


464,492 


436,464 


406,432 


445,474 


nig 


374 


392 


366 


339 


375 




-244 


-259 


-253 


-233 


-375 
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Table 3 



Scenario 1 










172(GeV) 


182 


194 


1 


0.117 pb 


0.155 


0.186 




(0.152) pb 


(0.198) 


(0.235) 


2 


0.132 


0.169 


0.198 




(0.169) 


(0.214) 


(0.248) 


3 


0.246 


0.267 


0.279 




(0.322) 


(0.344) 


(0.355) 


4 


0.117 


0.155 


0.186 




(0.152) 


(0.198) 


(0.235) 


5 


0.250 pb 


0.270 


0.283 




(0.322) pb 


(0.344) 


(0.356) 


eeWW 


0.177 


0.212 


0.225 
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Table 4 







-0.9 


-0.7 


-0.5 


-0.3 


-0.1 


0.1 


0.3 


0.5 


0.7 


0.9 


-0.9 


0.42 
(0.37) 


0.51 
(0.57) 


0.62 
(0.88) 


0.71 
(1.28) 


0.80 
(1.73) 


0.91 
(2.32) 


0.98 
(2.94) 


1.04 
(3.45) 


1.09 
(3.58) 


1.06 
(2.67) 


-0.7 


0.51 
(0.30) 


0.63 
(0.45) 


0.75 
(0.69) 


0.86 
(1.00) 


0.96 
(1.40) 


1.08 
(1.92) 


1.16 
(2.47) 


1.19 
(3.08) 


1.18 
(3.51) 




-0.5 


0.62 
(0.26) 


0.75 
(0.38) 


0.88 
(0.56) 


1.02 
(0.82) 


1.11 
(1.16) 


1.20 
(1.57) 


1.23 
(2.05) 


1.27 
(2.55) 






-0.3 


0.71 

(0.21) 


0.86 
(0.32) 


1.02 
(0.47) 


1.13 
(0.68) 


1.23 

(0.92) 


1.29 

(1.26) 


1.30 
(1.61) 








-0.1 


0.80 
(0.18) 


0.96 
(0.27) 


1.11 
(0.42) 


1.23 
(0.55) 


1.30 
(0.76) 


1.31 
(0.99) 










0.1 


0.91 
(0.16) 


1.08 
(0.23) 


1.20 
(0.38) 


1.29 
(0.45) 


1.31 
(0.58) 












0.3 


0.98 
(0.15) 


1.16 
(0.21) 


1.28 
(0.29) 


1.30 
(0.35) 














0.5 


1.04 
(0.15) 


1.19 
(0.19) 


1.27 
(0.23) 
















0.7 


1.09 
(0.15) 


1.18 
(0.17) 


















0.9 


1.06 
(0.15) 
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Table 5 



Scenarios 


172(GeV) 


182 


194 


6 


0.293 pb 


0.307 


0.312 


7 


0.132 


0.169 


0.198 


8 


7.22x10-2 


0.114 


0.151 


9 


8.54x10-2 


0.127 


0.162 


10 




4.39x10-3 


4.04x10-2 
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